Application of an integrated catchment-lake model approach for simulating effects of climate change on lake inputs and biogeochemistry.

Climate change is simultaneously affecting lakes and their catchments, resulting in altered runoff patterns in the catchment and modified mixing and biogeochemical dynamics in lakes. The effects of climate change in a catchment will eventually have an impact on the dynamics of a downstream water body as well. An integrated model would allow considering how changes in the watershed affect the lake, but coupled modelling studies are rare. In this study we integrate a catchment model (SWAT+) and a lake model (GOTM-WET) to obtain holistic predictions for Lake Erken, Sweden. Using five different global climate models, projections of climate, catchment loads and lake water quality for the mid and end of the 21st century have been obtained under two future scenarios (SSP 2-45 and SSP 5-85). Temperature, precipitation and evapotranspiration will increase in the future, overall resulting in an increase in water inflow to the lake. An increasing importance of surface runoff will also have consequences on the catchment soil, hydrologic flow paths, and the input of nutrients to the lake. In the lake, water temperatures will rise, leading to increased stratification and a drop in oxygen levels. Nitrate levels are predicted to remain unchanged, while phosphate and ammonium levels increase. A coupled catchment-lake configuration such as that illustrated here allows prediction of future biogeochemical conditions of a lake, including linking land use changes to changing lake conditions, as well as eutrophication and browning studies. Since climate affects both the lake and the catchment, simulations of climate change should ideally take into account both systems.

Assessing the effectiveness of potential best management practices for science-informed decision support at the watershed scale: The case of the Mar Menor coastal lagoon, Spain.

Coastal lagoons are ecosystems of high environmental importance but are quite vulnerable to human activities. The continuous inflow of pollutant loads can trigger negative impacts on the ecological status of these water bodies, which is contrary to the European Green Deal. One example is the Mar Menor coastal lagoon in Spain, which has experienced significant environmental degradation in recent years due to excessive external nutrient input, especially from non-point source (NPS) pollution. Mar Menor is one of the largest coastal lagoons of the Mediterranean region and a site of great ecological and socio-economic value. In this study, the highly anthropogenic and complex watershed of Mar Menor, known as Campo de Cartagena (1244 km2), was modelled with the Soil and Water Assessment Tool (SWAT) to analyse potential options for recovery of this unique system. The model was used to simulate several best management practices (BMP) proposed by recent Mar Menor regulations, such as vegetative filter strips, shoreline buffers, contour farming, removal of illegal agriculture, crop rotation management, waterway vegetation restoration, fertiliser management and greenhouse rainwater harvesting. Sixteen scenarios of individual and combined BMPs were analysed in this study. We found that, as individual measures, vegetative filter strips and contour farming were most effective in nutrient reduction: approximately 30 % for total nitrogen (TN) and 40 % for total phosphorus (TP). Moreover, waterway vegetation restoration showed the highest sediment (S) reduction at approximately 20 %. However, the combination of BMPs demonstrated clear synergistic effects, reducing S export by 38 %, TN by 67 %, and TP by 75 %. Selecting the most appropriate BMPs to be implemented at a watershed scale requires a holistic approach considering effectiveness in reducing NPS pollution loads and BMP implementation costs. Thus, we have demonstrated a way forward for enabling science-informed decision-making when choosing strategies to control NPS contamination at the watershed scale.

Impacts of land use, climate change and hydrological model structure on nitrate fluxes: Magnitudes and uncertainties.

Nitrate pollution and eutrophication are of increasing concern in agriculturally dominated regions, and with projected future climate changes, these issues are expected to worsen for both surface and groundwater. Changes in land use and management have the potential to mitigate some of these concerns. However, to what extent these changes will interact is unknown, and are associated with significant uncertainty. Here, we estimate nitrate fluxes and contributions of major uncertainty sources (variance decomposition analysis) affecting nitrate leaching from the root zone and river load from groundwater sources for an agricultural catchment in Denmark under future changes (2080-2099) in climate (four climate models) and land use (four land use scenarios). To investigate the uncertainty from impact model choice, two different agro-hydrological models (SWAT and DAISY-MIKE SHE) both traditionally used for nitrate impact assessments are used for projecting these effects. On average, nitrate leaching from the root zone increased by 55%-123% due to different climate models, while the impact of land use scenarios showed changes between -9% and 88%, with similar projections for river loads, while the worst-case combination of the three factors yielded a fivefold increase in nitrate transport. Thus, in the future, major land use changes will be necessary to mitigate nitrate pollution likely in combination with other measures such as advanced management and farming technologies and differentiated regulation. The two agro-hydrological models showed substantially different reaction patterns and magnitude of nitrate fluxes, and while the largest uncertainty source was the land use scenarios for both models, DAISY-MIKE SHE was to a higher degree affected by climate model choice. The dominating uncertainty source was found to be the agro-hydrological model; however, both uncertainties related to land use scenario and climate model were important, thus highlighting the need to include all influential factors in future nitrate flux impact studies.

Identifying key drivers of harmful algal blooms in a tributary of the Three Gorges Reservoir between different seasons: Causality based on data-driven methods.

Intense harmful algal blooms (HABs) can occur in the backwaters of tributaries supplying large-scale reservoirs. Due to the characteristics of process-based models and difficulties in modelling complex nonlinear processes, traditional models have difficulties disentangling the driving factors of HABs. In this study, we used data-driven methods (i.e., correlation analysis and machine-learning models) to identify the most important drivers of HABs in the Xiangxi River, a tributary of the Three Gorges Reservoir, China (2017-2018), for the dry season (from October to mid-April) and wet season (from April to September). We utilized the maximal information coefficient (MIC) combined with a time lag strategy and prior knowledge to quantitatively identify the driving variables of HABs. An extra trees regression (ETR) model was developed to assess the relative importance of causal variables driving algal blooms for the different periods. The results showed that water temperature was the most important driver for the duration of the study, followed by total nitrogen. Nitrogen had a stronger effect on algal blooms than phosphorus during both the wet and dry seasons. HABs were mainly affected by ammonia nitrogen in the wet season and by other forms of nitrogen in the dry season. In contrast, rather than the water temperature and nutrients, the operation of the Three Gorges Dam (difference between inflow and outflow discharge rate) was the most significant factor for algal blooms during the dry season, but its influence sharply declined during the wet season. This study showed that the key drivers of HABs can differ between seasons and suggests that HAB management should take seasonality into account.

Subsurface hydrological processes and groundwater residence time in a coastal alluvium aquifer: Evidence from environmental tracers (δ18O, δ2H, CFCs, 3H) combined with hydrochemistry.

As an important part of the water cycle, the hydrologic process and chemical compositions of groundwater have changed significantly due to the joint influence of climate change and human activities. Groundwater salinization becomes a serious threat to water security in coastal areas. In order to assess the relationships between surface water, groundwater and seawater in the coastal plain, we performed a synthesis study based on hydrochemical-isotopic data, hydro-dynamical records and environmental tracers. Deuterium and oxygen isotopes and water chemical indicators were used to identify pollution status, salt sources and migration processes. Radioactive isotopes and gaseous tracers were used to obtain reasonable groundwater age. With the help of multi-tracer approach, the surface-groundwater interaction, salinization of groundwater and nitrate pollution were identified in the Yang-Dai River plain, northern China. The estimated groundwater ages determined from chlorofluorocarbons (CFCs) and tritium (3H) ranges from 18 to 41 years in this area, suggesting a modern groundwater circulation. The spatial distribution of the groundwater age varies significantly due to horizontal hydrogeological heterogeneity. The total dissolved solids (TDS) content of the groundwater near the Well Field (average: 970 mg/L) was higher than the TDS values in samples derived from places located at an equivalent distance to the coastal line (average is 708 mg/L), which resulted from the vertical seawater intrusion through river channels and pollutant inputs from agriculture activities. The nitrate concentrations in groundwater were elevated up to 271 mg/L and increased with increasing groundwater age, which was another water environment problem that should be solved urgently but lacks sufficient attention for years. This study provides a conceptual model with a number of comparable hydrochemical information, which is significant for regional pollution control and water resources management.

Quantifying the effects of climate change on hydrological regime and stream biota in a groundwater-dominated catchment: A modelling approach combining SWAT-MODFLOW with flow-biota empirical models.

Climate change may affect stream ecosystems through flow regime alterations, which can be particularly complex in streams with a significant groundwater contribution. To quantify the impacts of climate change on hydrological regime and subsequently the stream biota, we linked SWAT-MODFLOW (A model coupling the Soil and Water Assessment Tool and the Modular Finite-difference Flow Model) with flow-biota empirical models that included indices for three key biological taxonomic identities (fish, macroinvertebrates and macrophytes) and applied the model-complex to a groundwater-dominated catchment in Denmark. Effects of predicted climate change towards the end of this century relative to the reference period (1996-2005) were tested with two contrasting climate change scenarios of different greenhouse gas emissions (Representative Concentration Pathway 2.6 (RCP 2.6) and RCP 8.5) and analysed for all subbasins grouped into streams of three size classes. The total water yield in the catchment did not change significantly (-1 ± 4 (SD) mm yr-1) from the baseline in the RCP2.6 scenario, while it increased by 9 ± 11 mm yr-1 in the RCP8.5 scenario. The three stream size classes underwent different alterations in flow regime and also demonstrated different biotic responses to climate change. All large and some small streams were impacted most heavily by the climate change, where fish and macrophyte indices decreased up to 14.4% and 11.2%, respectively, whereas these indices increased by up to 14.4% and 6.0%, respectively, in the medium and some small streams. The climate change effects were, as expected, larger in the RCP8.5 scenario than in the RCP2.6 scenario. Our study is the first to quantify the impacts of streamflow alterations induced by climate change on stream biota beyond specific species.

Predicting ecosystem state changes in shallow lakes using an aquatic ecosystem model: Lake Hinge, Denmark, an example.

In recent years, considerable efforts have been made to restore turbid, phytoplankton-dominated shallow lakes to a clear-water state with high coverage of submerged macrophytes. Various dynamic lake models with simplified physical representations of vertical gradients, such as PCLake, have been used to predict external nutrient load thresholds for such nonlinear regime shifts. However, recent observational studies have questioned the concept of regime shifts by emphasizing that gradual changes are more common than sudden shifts. We investigated if regime shifts would be more gradual if the models account for depth-dependent heterogeneity of the system by including the possibility of vertical gradients in the water column and sediment layers for the entire depth. Hence, bifurcation analysis was undertaken using the 1D hydrodynamic model GOTM, accounting for vertical gradients, coupled to the aquatic ecosystem model PCLake, which is implemented in the framework for aquatic biogeochemical modeling (FABM). First, the model was calibrated and validated against a comprehensive data set covering two consecutive 7-yr periods from Lake Hinge, a shallow, eutrophic Danish lake. The autocalibration program Auto-Calibration Python (ACPy) was applied to achieve a more comprehensive adjustment of model parameters. The model simulations showed excellent agreement with observed data for water temperature, total nitrogen, and nitrate and good agreement for ammonium, total phosphorus, phosphate, and chlorophyll a concentrations. Zooplankton and macrophyte coverage were adequately simulated for the purpose of this study, and in general the GOTM-FABM-PCLake model simulations performed well compared with other model studies. In contrast to previous model studies ignoring depth heterogeneity, our bifurcation analysis revealed that the spatial extent and depth limitation of macrophytes as well as phytoplankton chlorophyll-a responded more gradually over time to a reduction in the external phosphorus load, albeit some hysteresis effects still appeared. In a management perspective, our study emphasizes the need to include depth heterogeneity in the model structure to more correctly determine at which external nutrient load a given lake changes ecosystem state to a clear-water condition.

Assessing the impacts of groundwater abstractions on flow regime and stream biota: Combining SWAT-MODFLOW with flow-biota empirical models.

Assessing the impacts of groundwater abstractions on stream ecosystems is crucial for developing water planning and regulations in lowland areas that are highly dependent on groundwater, such as Denmark. To assess the effects of groundwater abstractions on flow regime and stream biota in a lowland groundwater-dominant catchment, we combined the SWAT-MODFLOW model with flow-biota empirical models including indices for three key biological taxonomic identities (fish, macroinvertebrates, and macrophytes). We assessed the effects of the current level of abstractions and also ran a scenario for assessing the effect of extreme groundwater abstractions (pumping rates of the drinking water wells were increased by 20 times in one subbasin of the catchment). Three subbasin outlets representing stream segments of different sizes were used for this evaluation. Current groundwater abstraction level had only minor impacts on the flow regime and stream biotic indices at the three subbasin outlets. The extreme abstractions, however, led to significant impacts on the small stream but had comparatively minor effects on the larger streams. The fish index responded most negatively to the groundwater abstractions, followed by the macrophyte index, decreasing, respectively, by 23.5% and 11.2% in the small stream in the extreme groundwater abstraction scenario. No apparent impact was found on macroinvertebrates in any of the three subbasin outlets. We conclude that this novel approach of a combined modelling system is a useful tool to quantitatively assess the effects of groundwater abstractions on stream biota and thereby support water planning and regulations related to groundwater abstractions. We highlight the need for developing improved biotic models that target specifically small headwater streams, which are often most affected by water abstraction.

Effects of changes in land use and climate on aquatic ecosystems: Coupling of models and decomposition of uncertainties.

To analyse the potential future ecological state of estuaries located in the temperate climate (here exemplified with the Odense Fjord estuary, Denmark), we combined end-of-the-century climate change projections from four different climate models, four contrasting land use scenarios ("Agriculture for nature", "Extensive agriculture", "High-tech agriculture" and "Market driven agriculture") and two different eco-hydrological models. By decomposing the variance of the model-simulated output from all scenario and model combinations, we identified the key sources of uncertainties of these future projections. There was generally a decline in the ecological state of the estuary in scenarios with a warmer climate. Strikingly, even the most nature-friendly land use scenario, where a proportion of the intensive agricultural area was converted to forest, may not be enough to counteract the negative effects of a future warmer climate on the ecological state of the estuary. The different land use scenarios were the most significant sources of uncertainty in the projections of future ecological state, followed, in order, by eco-hydrological models and climate models, albeit all three sources caused high variability in the simulated outputs. Therefore, when projecting the future state of aquatic ecosystems in a global warming context, one should at the very least consider to evaluate an ensemble of land use scenarios (nutrient loads) but ideally also include multiple eco-hydrological models and climate change projections. Our study may set precedence for future attempts to predict and quantify uncertainties of model and model input ensembles, as this will likely be key elements in future tools for decision-making processes.

Modelling the fate and transport of Cryptosporidium, a zoonotic and waterborne pathogen, in the Daning River watershed of the Three Gorges Reservoir Region, China.

Oospores of Cryptosporidium, a waterborne pathogen of great concern, are widely distributed in surface waters in China and pose a threat to human health. This study seeks to explore the spatio-temporal patterns of Cryptosporidium concentrations. We focus on the Daning River watershed (4166 km2) of the Three Gorges Reservoir Region (TGR) during the period 2008 to 2013 and use the SWAT (Soil and Water Assessment Tool) model to test two mitigation scenarios. Based on data on animal husbandry, population, agriculture and WWTPs (wastewater treatment plants), Cryptosporidium transport in the Daning River watershed was simulated using a calibrated hydrological and sediment transport model. Our model results showed that the average annual concentration of oocysts in the whole watershed was 9.5 oocysts/10L, but high spatial variability occurred, ranging from 0.7 to 33.4 oocysts/10L. Highest monthly mean oocysts concentrations at the outlets of the sub-basins were found at high runoff and high fertilization or at the lowest flow, while minimum monthly mean oocysts concentrations were recorded at high runoff only. A model parameter sensitivity analysis showed that the Cryptosporidium soil partitioning coefficient (BACTKDQ) and the temperature adjustment factor for Cryptosporidium die-off (THBACT) were the only two sensitive parameters among the microbial parameters. The construction of multiple WWTPs throughout the watershed and composting of 50% of the feces from rural citizens and livestock up to 56 days before its application as fertilizer could significantly reduce the concentration of oocysts. Our Cryptosporidium transport model and simulation results may assist in the establishment of better pollution control countermeasures in the Daning River and other similar watersheds.

Modeling the effects of climatic and land use changes on phytoplankton and water quality of the largest Turkish freshwater lake: Lake Beysehir.

Climate change and intense land use practices are the main threats to ecosystem structure and services of Mediterranean lakes. Therefore, it is essential to predict the future changes and develop mitigation measures to combat such pressures. In this study, Lake Beysehir, the largest freshwater lake in the Mediterranean basin, was selected to study the impacts of climate change and various land use scenarios on the ecosystem dynamics of Mediterranean freshwater ecosystems and the services that they provide. For this purpose, we linked catchment model outputs to the two different processed-based lake models: PCLake and GLM-AED, and tested the scenarios of five General Circulation Models, two Representation Concentration Pathways and three different land use scenarios, which enable us to consider the various sources of uncertainty. Climate change and land use scenarios generally predicted strong future decreases in hydraulic and nutrient loads from the catchment to the lake. These changes in loads translated into alterations in water level as well as minor changes in chlorophyll a (Chl-a) concentrations. We also observed an increased abundance of cyanobacteria in both lake models. Total phosphorus, temperature and hydraulic loading were found to be the most important variables determining cyanobacteria biomass. As the future scenarios revealed only minor changes in Chl-a due to the significant decrease in nutrient loads, our results highlight that reduced nutrient loading in a warming world may play a crucial role in offsetting the effects of temperature on phytoplankton growth. However, our results also showed increased abundance of cyanobacteria in the future may threaten ecosystem integrity and may limit drinking water ecosystem services. In addition, extended periods of decreased hydraulic loads from the catchment and increased evaporation may lead to water level reductions and may diminish the ecosystem services of the lake as a water supply for irrigation and drinking water.

Quantifying the combined effects of land use and climate changes on stream flow and nutrient loads: A modelling approach in the Odense Fjord catchment (Denmark).

Water pollution and water scarcity are among the main environmental challenges faced by the European Union, and multiple stressors compromise the integrity of water resources and ecosystems. Particularly in lowland areas of northern Europe, high population density, flood protection and, especially, intensive agriculture, are important drivers of water quality degradation. In addition, future climate and land use changes may interact, with uncertain consequences for water resources. Modelling approaches have become essential to address water issues and to evaluate ecosystem management. In this work, three multi-stressor future storylines combining climatic and socio-economic changes, defined at European level, have been downscaled for the Odense Fjord catchment (Denmark), giving three scenarios: High-Tech agriculture (HT), Agriculture for Nature (AN) and Market-Driven agriculture (MD). The impacts of these scenarios on water discharge and inorganic and organic nutrient loads to the streams have been simulated using the Soil and Water Assessment Tool (SWAT). The results revealed that the scenario-specific climate inputs were most important when simulating hydrology, increasing river discharge in the HT and MD scenarios (which followed the high emission 8.5 representative concentration pathway, RCP), while remaining stable in the AN scenario (RCP 4.5). Moreover, discharge was the main driver of changes in organic nutrients and inorganic phosphorus loads that consequently increased in a high emission scenario. Nevertheless, both land use (via inputs of fertilizer) and climate changes affected the nitrate transport. Different levels of fertilization yielded a decrease in the nitrate load in AN and an increase in MD. In HT, however, nitrate losses remained stable because the fertilization decrease was counteracted by a flow increase. Thus, our results suggest that N loads will ultimately depend on future land use and management in an interaction with climate changes, and this knowledge is of utmost importance for the achievement of European environmental policy goals.

Future water availability in the largest freshwater Mediterranean lake is at great risk as evidenced from simulations with the SWAT model.

Inter- and intra-annual water level fluctuations and changes in water flow regime are intrinsic characteristics of Mediterranean lakes. Additionally, considering climate change projections for the water-limited Mediterranean region, increased air temperatures and decreased precipitation are anticipated, leading to dramatic declines in lake water levels as well as severe water scarcity problems. The study site, Lake Beysehir, the largest freshwater lake in the Mediterranean basin, is - like other Mediterranean lakes - threatened by climatic changes and over-abstraction of water for irrigated crop farming. Therefore, implementation of strict water level management policies is required. In this study, an integrated modeling approach was used to predict the future water levels of Lake Beysehir in response to potential future changes in climate and land use. Water level estimation was performed by linking the catchment model Soil and Water Assessment Tool (SWAT) with a Support Vector Regression model (ε-SVR). The projected increase in temperature and decrease in precipitation based on the climate change models led to an enhanced potential evapotranspiration and reduced total runoff. On the other hand, the effects of various land use scenarios within the catchment appeared to be comparatively insignificant. According to the ε-SVR model results, changes in hydrological processes caused a water level reduction for all scenarios. Moreover, the MPI-ESM-MR General Circulation Model outputs produced the most dramatic results by predicting that Lake Beysehir may dry out by the 2040s with the current outflow regime. The results indicate that shallow Mediterranean lakes may face a severe risk of drying out and losing their ecosystem values in the near future if the current intensity of water abstraction is not reduced. In addition, the results also demonstrate that outflow management and sustainable use of water sources are vital to sustain lake ecosystems in water-limited regions.

Application of a three-dimensional water quality model as a decision support tool for the management of land-use changes in the catchment of an oligotrophic lake.

While expansion of agricultural land area and intensification of agricultural practices through irrigation and fertilizer use can bring many benefits to communities, intensifying land use also causes more contaminants, such as nutrients and pesticides, to enter rivers, lakes, and groundwater. For lakes such as Benmore in the Waitaki catchment, South Island, New Zealand, an area which is currently undergoing agricultural intensification, this could potentially lead to marked degradation of water clarity as well as effects on ecological, recreational, commercial, and tourism values. We undertook a modeling study to demonstrate science-based options for consideration of agricultural intensification in the catchment of Lake Benmore. Based on model simulations of a range of potential future nutrient loadings, it is clear that different areas within Lake Benmore may respond differently to increased nutrient loadings. A western arm (Ahuriri) could be most severely affected by land-use changes and associated increases in nutrient loadings. Lake-wide annual averages of an eutrophication indicator, the trophic level index (TLI) were derived from simulated chlorophyll a, total nitrogen, and total phosphorus concentrations. Results suggest that the lake will shift from oligotrophic (TLI = 2-3) to eutrophic (TLI = 4-5) as external loadings are increased eightfold over current baseline loads, corresponding to the potential land-use intensification in the catchment. This study provides a basis for use of model results in a decision-making process by outlining the environmental consequences of a series of land-use management options, and quantifying nutrient load limits needed to achieve defined trophic state objectives.

Effects of climate and nutrient load on the water quality of shallow lakes assessed through ensemble runs by PCLake.

Complex ecological models are used to predict the consequences of anticipated future changes in climate and nutrient loading for lake water quality. These models may, however, suffer from nonuniqueness in that various sets of model parameter values may yield equally satisfactory representations of the system being modeled, but when applied in future scenarios these sets of values may divert considerably in their simulated outcomes. Compilation of an ensemble of model runs allows us to account for simulation variability arising from model parameter estimates. Thus, we propose a new approach for aquatic ecological models creating a more robust prediction of future water quality. We used our ensemble approach in an application of the widely used PCLake model for Danish shallow Lake Arreskov, which during the past two decades has demonstrated frequent shifts between turbid and clear water states. Despite marked variability, the span of our ensemble runs encapsulated 70­90% of the observed variation in lake water quality. The model exercise demonstrates that future warming and increased nutrient loading lead to lower probability of a clear water, vegetation-rich state and greater likelihood of cyanobacteria dominance. In a 6.0°C warming scenario, for instance, the current nutrient loading of nitrogen and phosphorus must be reduced by about 75% to maintain the present ecological state of Lake Arreskov, but even in a near-future 2.0°C warming scenario, a higher probability of a turbid, cyanobacteria-dominated state is predicted. As managers may wish to determine the probability of achieving a certain ecological state, our proposed ensemble approach facilitates new ways of communicating future stressor impacts.

Watershed land use effects on lake water quality in Denmark.

Mitigating nutrient losses from anthropogenic nonpoint sources is today of particular importance for improving the water quality of numerous freshwater lakes worldwide. Several empirical relationships between land use and in-lake water quality variables have been developed, but they are often weak, which can in part be attributed to lack of detailed information about land use activities or point sources. We examined a comprehensive data set comprising land use data, point-source information, and in-lake water quality for 414 Danish lakes. By excluding point-source-influenced lakes (n = 210), the strength in relationship (R2) between in-lake total nitrogen (TN) and total phosphorus (TP) concentrations and the proportion of agricultural land use in the watershed increased markedly, from 10-12% to 39-42% for deep lakes and from 10-12% to 21-23% for shallow lakes, with the highest increase for TN. Relationships between TP and agricultural land use were even stronger for lakes with rivers in their watershed (55%) compared to lakes without (28%), indicating that rivers mediate a stronger linkage between landscape activity and lake water quality by providing a "delivery" mechanism for excess nutrients in the watershed. When examining the effect of different near-freshwater land zones in contrast to the entire watershed, relationships generally improved with size of zone (25, 50, 100, 200, and 400 m from the edge of lake and streams) but were by far strongest using the entire watershed. The proportion of agricultural land use in the entire watershed was best in explaining lake water quality, both relative to estimated nutrient surplus at agricultural field level and near-lake land use, which somewhat contrasts typical strategies of management policies that mainly target agricultural nutrient applications and implementation of near-water buffer zones. This study suggests that transport mechanisms within the whole catchment are important for the nutrient export to lakes. Hence, the whole watershed should be considered when managing nutrient loadings to lakes, and future policies should ideally target measures that reduce the proportion of cultivated land in the watershed to successfully improve lake water quality.